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Characterization of real- world complex systems increasingly involves the study of their topological structure 
using graph theory. Among global network properties, small-world property, consisting in existence of rel- 
atively short paths together with high clustering of the network, is one of the most discussed and studied. 
When dealing with coupled dynamical systems, links among units of the system are commonly quantified by 
a measure of pairwise statistical dependence of observed time series (functional connectivity) . We argue that 
the functional connectivity approach leads to upwardly biased estimates of small-world characteristics (with 
respect to commonly used random graph models) due to partial transitivity of the accepted functional connec- 
tivity measures such as the correlation coefficient. In particular, this may lead to observation of small-world 
characteristics in connectivity graphs estimated from generic randomly connected dynamical systems. The 
ubiquity and robustness of the phenomenon is documented by an extensive parameter study of its manifes- 
tation in a multivariate linear autoregressive process, with discussion of the potential relevance for nonlinear 
processes and measures. 

PACS numbers: 89.75.Fb,07.05.Kf,05.45.Tp,87.19.L- 
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In the field of complex systems study, new mea- 
surement and computational resources have lead 
to increased interest in analysis of large networks. 
These networks are observed across many disci- 
plines spanning from social sciences through bi- 
ology to climate research. For characterization of 
the structure of these networks, graph-theoretical 
measures have proven to be useful. These charac- 
teristics capture some global features of the net- 
work topology such as the density or level of clus- 
tering (tendency of neighbors of a node to be also 
neighbors to each other) as well as specific roles of 
important nodes serving as 'hubs' in the network. 
Among the interesting properties of many real- 
world networks belongs the small-world property, 
a global property of a network characterized by 
a relatively high level of clustering while con- 
serving on average short paths among nodes of 
the network, compared to a random network of 
corresponding density. This small-world prop- 
erty has been related to some convenient prop- 
erties of the network including efficiency of infor- 
mation transfer, and therefore reports of small- 
world architecture in real-world networks have 
received much attention. In some complex dy- 
namical systems, including global climate or hu- 
man brain, the knowledge of physical connections 
among its subsystems is far from perfect, and 
therefore other methods of characterizing interac- 
tions among these have been extensively applied. 
In particular, interactions among the areas are 
commonly quantified by a dependence measure 
such as the linear correlation between the local 
time series of variables of interest. This gives 
rise to the so-called functional connectivity ma- 



trix of a system. Applying the graph-theoretical 
approach to functional connectivity matrices has 
lead to reports of small-world properties of many 
real-world systems. However, as we document in 
this report, even for a simplistic dynamical sys- 
tem with linear dynamics and random coupling 
matrix the functional connectivity approach gen- 
erates networks with small-world characteristics. 
These spurious detections of small-world topol- 
ogy is related to partial transitivity of functional 
connectivity measures such as the correlation co- 
efficient. 



I. INTRODUCTION 

Characterization of complex systems commonly in- 
cludes the study of their structure using graph theory. 
This typically involves identification of the systems sub- 
units (nodes of a network) and assessment of existence 
(or strength) of pair- wise relations among those, leading 
to representation of the system by a (weighted) graph. 
The local, mesoscale or global topology of the graph (or 
alterations thereof) are subsequently studied using vari- 
ous graph-theoretical measures with a goal of identifying 
systems properties key to its function. Among the most 
widely discussed properties are modularity, scale-freeness 
or small- world topology. 

Since the paradigmatic publication of Watts and Stro- 
gatz in Nature 1 , the small- world property entailing rela- 
tively short graph paths and high clustering has received 
much attention in many application areas dealing with 
complex systems. This includes such diverse fields as 
ncuroscience 2 and climate science 3 . 
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When working with such complex dynamical systems, 
for both practical and theoretical reasons, the links 
among nodes are commonly quantified by the dependence 
of the observed time series rather than the underlying 
physical or coupling network of connections. In the neu- 
roscience community, this corresponds to the distinction 
between functional and structural (anatomical) connec- 
tivity. We adopt this useful terminology in this paper and 
suggest the general term functional connectivity approach 
to construction of graph representation of complex dy- 
namical system by deriving the links through quantifica- 
tion (and potential thresholding) of the pair-wise statis- 
tical dependence measure (typically correlation) among 
the local time series of activity or other key variable of 
the nodes of the studied system. 

As the exact values of quantitative graph-theoretical 
indices crucially depend on parameters such as network 
size or density, the interpretation of graph-theoretical 
properties of real-world networks is usually based on the 
comparison with 'corresponding' random graph models 
of Erdos-Renyi 4 or Maslov-Sneppen 5 . 

Some complications with the outlined analysis ap- 
proach have recently been suggested in literature, partic- 
ularly discussing the biasing effects related to sampling 
problems such as spatial oversampling 6 and finite size 
temporal samples and their autocorrelation 7 . 

However, we argue that there is a more fundamen- 
tal problem with intepretation of the graph-theoretical 
properties of functional connectivity matrices. In partic- 
ular, functional connectivity matrices are biased towards 
a specific structure due to their construction method; ir- 
respective of particular time series length and/or sam- 
pling parameters. We focus particularly on functional 
connectivity matrices constructed using linear correla- 
tion as a measure of dependence. A specific example 
of an in-built bias of a linear correlation matrix is its 
'weak' transitivity property: for any three random vari- 
ables X, Y, Z a strong positive correlation between two 
pairs of them implies a positive correlation within the 
third pair, including specific 'hard limitations', such as 
that p 2 XY + Pyz > 1 implies pxz > 0, i.e. positivity of 
the third correlation coefficient (for a proof of a general 
form of this inequality see 8 ). Such implicit dependence 
among the entries of the correlation matrix has been com- 
monly overlooked in the interpretation of graph analysis 
of functional connectivity matrices. 

In the presented study, we focus on the effects of this 
dependence on the graph-theoretical properties of the 
functional connectivity matrices. In particular, we show 
that there are tendencies towards specific graph struc- 
tures in functional connectivity matrices computed from 
activity time-series of randomly structurally connected 
networks. A prime example of such an effect is the in- 
creased level of clustering, further affecting the estimates 
of small- world indices, potentially leading to serious mis- 
interpretation of real- world data. 



II. METHODS 

A. Functional connectivity 

Functional connectivity is most commonly quantified 
by the linear (Pearson's) correlation coefficient. While 
also non-linear connectivity measures are in use, for data 
with approximately Gaussian distribution is linear cor- 
relation practically sufficient 9 ' 10 . The correlation ma- 
trix can be subsequently transformed into an unweighted 
graph ('binarized') by choosing a threshold and assign- 
ing links only to pairs of nodes with over-threshold cor- 
relation. In practice, the threshold is commonly chosen 
adaptively to achieve a predefined density of the result- 
ing graph, although statistical testing may be applied as 
well 11 - 13 . 



B. Graph-theoretical characteristics 

Formally, in the graph-theoretical approach a network 
is represented by an unweighted graph G = (V, E), where 
V is the set of nodes of G, N = #V is the number of 
nodes and E C (^) is the set of the edges of G. In some 
cases a separate function w : E —¥ R can be used to 
extend the graph concept by defining a weighted graph 
G = (V, E, w), however throughout the current paper we 
only deal with unweighted graphs. 

For any pair of nodes i, j ^ V we define an edge indi- 
cator aij = 1 if and only if {i,j} £ E (the graph con- 
tains an edge linking nodes i and j), otherwise Ojj = 0, 
node degree as ki = X^gy a i,j-> an d finally a path of 
length t from node i to node j as a sequence of edges 
{to = i,mi}, {mi,m 2 }, . . ., {m<_i,m< = j} where all 
TOj G V are distinct nodes of the graph. The length of 
the shortest path between nodes i and j is denoted as 
dij . 

A graph G can be characterized by its global proper- 
ties, including the average path length 

and the clustering coefficient 

iev 

where a denotes the local clustering coefficient defined 
for ki > 2 as 

Cl ~ ki(ki-l) [6) 

and for ki £ {0, 1} as Cj = 0. 

Note that for disconnected graphs, some pairs of nodes 
are not connected by a path and the respective dij is not 
well defined (or set to infinity). In this case we compute 
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the average path length by averaging only over pairs for 
which a finite path exists, see also the discussion section. 

In their original paper, Watts and Strogatz 1 suggested 
the term "small-world" for networks that have a similar 
average path length, but an increased clustering coeffi- 
cient compared to a corresponding random graph (which 
corresponds to a relative average path length and clus- 
tering coefficient A = L L - > 1, 7 = c c - ^> 1) and 
also provided a simple generative modef for such net- 
works. The small-world property has been more recently 
proposed to be summarized in the small-world index 

Although generally accepted, the interpretation of the 
small-world index has been critically discussed for var- 
ious reasons including the potential confounding effects 
of spatial oversampling of the underlying system 6 , finite- 
time estimation effects 7 or unclear interpretation of paths 
in the functional connectivity graphs 15 . In the following, 
we show that the functional connectivity approach with 
a common choice of dependence measure leads to a sub- 
stantial upward bias in the small-world index a. Under- 
standing this phenomenon is key for accurate interpreta- 
tion of experimental findings and their comparison across 
methodologies. 

For simplicity, consider an autoregressive process of 
order 1 (AR(1)): 

X t = c + AX t _ 1 +e t , (4) 

where c is a N x I vector of constants, A is a, N x N 
matrix and e* is a N x I vector of error terms. For 
simplicity we choose c = Ojv.i and e t ~ Af(0,l) and 
A = s(SC + aI)/X max , with the symmetric binary struc- 
tural connectivity matrix SC = SC (N, p) generated as a 
realization of the Erdos-Renyi model G(N,p). In partic- 
ular, each nondiagonal entry of the N x N SC matrix is 
assigned randomly and independently either with value 
1 (edge exists, with probability p) or value (edge does 
not exist, with probability 1 — p); diagonal elements are 
set to 0. I here denotes the identity matrix. The pa- 
rameter s G (0,1) modulates the relative strength of the 
autoregressive and noise terms in (4) with the normaliza- 
tion by Xmax ~ the largest (in absolute value) eigenvalue 
of the matrix SC + al. The parameter a varies the re- 
lation between the autocorrelation and cross-correlation 
component of the AR process (4). 



C. Finite time series example 

As a motivational example, we generate a finite sample 
of stochastic process with length T = 300 with parameter 
setting psc = 0.1, iV = 100, s = 0.1, a = 2, and compute 
the corresponding functional connectivity matrix FC by 
binarizing the sample Pearson correlation matrix. The 
binarization threshold is chosen such that the density 
Pfc of the binarized functional connectivity matrix is 
equal to the density psc of the structural connectivity 



matrix (diagonal elements of the FC matrix are first set 
to zero). The respective matrices are shown in Figure 1. 

A visual inspection shows that while the entries of the 
structural connectivity matrix are random and mutually 
independent, the functional connectivity matrix shows a 
specific structure. This can be quantified by the graph- 
theoretical measures. In this particular realization, we 
have L s = 2.157, L F = 2.308, C s = 0.1081, C F = 0.2355. 
As SC and FC have the same densities, SC is effectively 
a realization of the Erdos-Renyi model corresponding to 
the density of FC, and we obtain the relative graph mea- 
sures: A = 1.07,7 = 2.18,(7 = 2.04. The values indicate 
increased clustering and approximately conserved aver- 
age path length with respect to a corresponding random 
graph. Together this signifies a small-world like struc- 
ture of the functional connectivity matrix, even though 
the coupling structure of the generating system is com- 
pletely random. 

In view of the results of this simple experiment, it is 
important to ask what the significance of the findings of 
increased clustering and small- world structure in depen- 
dence matrices of real- world data is. 



D. Asymptotic behavior of correlation matrix: Parametric 
study settings 

In the following we focus on investigating this effect in 
a more detail, studying its strength as a function of some 
of the most relevant parameters of the underlying pro- 
cess. For the sake of tractability we still limit ourselves 
to AR(1) process, but vary the following parameters: size 
of the network N, density psc of the coupling matrix SC, 
balance of the autocorrelation and cross-correlation (by 
parameter a), balance between the autoregressive and 
noise terms by parameter s and the level of threshold- 
ing by changing the required density pfc of the matrix 
FC. While the strength of the effect is also dependent on 
the length of the sample T, the finite size of the sample 
is not crucial, as the upward bias does not vanish even 

asymptotically. The covariance matrix S for the AR(1) 

00 

process (4) is given by the infinite sum E = ^ A 1 A 1 

i=0 

(for derivation see e.g. 16 ), which due to symmetry of the 
matrix A = A T leads to a Neumann series that con- 
verges to £ = (I — A 2 )^ 1 . From this expression the 
correlation matrix is easily obtained by trivial normal- 
izations using diagonal elements of S. Thus the choice of 
the rather trivial AR(1) model gives us the possibility to 
study asymptotic behavior of correlation matrices with- 
out additional effects of sampling or series length which 
are not negligible when the correlation matrices are esti- 
mated from time series. 

Notwithstanding the dimensionality of the investi- 
gated parameter space, the parameter ranges were 
chosen to cover reasonably the parameter space, in 
particular we used values N £ {50,200,500}, s <E 
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FIG. 1. An example of binary functional connectivity matrix (right) generated from random structural connectivity matrix 
(left) by thresholding the correlation matrix of AR-model generated time series (center, light shades of gray indicate higher 
correlation values). Network with N = 100 nodes shown. Note that the functional connectivity matrix shows a specific structure 
although the entries of the generating structural connectivity matrix were chosen randomly. See text for further details. 



{0.2,0.5,0.75,0.9,0.99}, a e {0,1}. We further varied 
Psc and Pfc logarithmically in 24 steps within the (0, 1) 
interval - more exactly both variables are defined as 2" 
where n is an arithmetic progression from to —6.9 with 
step —0.3. The lowest density was therefore smaller than 
0.01. 

For robustness of evidence, for each parameter setting 
we compute 20 independent realizations of the coupling 
matrix, and each of the resulting matrices FC is com- 
pared to its own corresponding realization of the random 
Erdos-Rcnyi matrix G (and secondarily also to Maslov- 
Sneppen random graph model 5 , see discussion). The 
computations were carried out using the NDW-Graph 
Toolbox (http://ndw.cs.cas.cz/software/ndw-graph), a 
C++/MATLAB toolbox for complex network analysis 
developed within the authors' NDW group and available 
to the public under GPL license. 



III. RESULTS 

The simulations have shown that the studied effect 
(small-world property of functional connectivity matrix 
of a randomly connected dynamical system) is present 
throughout the covered parameter space, although with 
a variable strength. A representative selection of results 
is shown in Figure 2, obtained for network size N = 500 
nodes, s = 0.75 and a — 1. The results for other investi- 
gated parameter choices are qualitatively similar and are 
summarized for completeness in Figures 1 and 2 of the 
Supplemental Material 17 . 

Within the studied parameter range, we have observed 
values of small-world index a up to the order of hun- 
dreds. The effect has also proven statistically robust with 
respect to different realizations of the structural matrix 
SC. In particular, we observed relatively low spread of 
the a values, with a > 1 in all 20 assessed realizations 
of the process (4) for overwhelming majority of param- 
eter vector values. This corresponds to robust statis- 



tical significance in most cases (p- values < 10 -5 , sign 
test of hypothesis of median equal to 1, no correction for 
multiple comparisons; similar results obtained for t-test). 
The only exceptions were observed for the case of exactly 
equal densities of structural and functional connectivity 
matrix, when this commong density was very low (e.g. 
only for pfc = Psc ^ 0.03 for the specific settings in 
Figure 2), where the a values were relatively close to 1; 
this special case is discussed later. 

In general, a increases with increasing thresholding 
(that is decreasing density pfc of the FC matrix). The 
dependence on the density of the underlying structural 
connectivity psc is not monotonous (see Figure 2). Ap- 
proximately, the effect is weakest if the density of the 
binarized functional connectivity matrix is the same as 
of the underlying structural connectivity matrix, but the 
depth of this minimum further depends on the value of 
parameter a. The values of a are comparable for both 
a = 1 and a = 0, with the exception of the above- 
described non-monotonicity being more pronounced with 
a = 1. With other parameters fixed, network size N and 
strength of coupling s do not affect a strongly within the 
sampled parameter region (see Supplemental Material 17 
Figure 1). 

However, for small networks (especially N = 50) and 
small densities of the FC matrix ppc < 0.1, the a is often 
not well defined, mostly due to zero values of the clus- 
tering coefficient of the corresponding random matrix. 
Moreover, in a part of the parameter space, the result- 
ing graphs are not connected. This generally happens 
for small density, however the exact threshold is further 
modulated by network size and other network parame- 
ters. 

The interpretation of some graph-theoretical measures 
for disconnected graphs and their comparison to con- 
nected graphs is not straightforward. In practice, for 
disconnected graphs the average path length is computed 
only for the largest component of the graph 1 , or rede- 
fined as the mean of the finite elements of the distance 
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FIG. 2. Average values of the small-world index a of the bi- 
narized functional connectivity matrix FC generated by AR 
model with random structural connectivity matrix SC. In 
most of the parameter space, a>l, suggesting small-world 
topology. Note marked dependence on values of the density 
of SC and FC (the latter modulated by adjusting the bina- 
rization threshold). Other parameter values: s = 0.75, a = 
1, iV = 500. Logarithmic scale used both for the ordinate 
and the abscissa. Undefined values (for SC = 1; FC / 1) are 
shown as (in white). (Results for other parameter settings 
are available in the Supplemental Material 17 .) 



matrix, omitting the unconnected node pairs. In the cur- 
rent study we adopt the latter approach, used previously 
in similar context 6 . 

However, from theoretical perspective it could be ar- 
gued that this could potentially have biasing effects on 
the small- world index estimates, supporting a conserva- 
tive approach of only dealing with connected graphs. To 
check that our results are not driven only by loss of con- 
nectedness of the graphs with decreasing density or the 
related redefinition of average path, we plot separately 
the results obtained from connected graphs. While this 
severely limits the parameter space, even within this area 
a values attained values up to ~ 10 (see Supplemental 
Material 17 , Figure 3). 

As the admissible area of the parameter space 
monotonously grows with the network size, this docu- 
ments that even the conservative connectedness require- 
ment in principle does not limit the values of a and there- 
fore the strength of the confounding effect. 

The contributions of relative clustering coefficient 7 
and relative average path length A to the observed small- 
world property is documented in Figures 3 and 4. The 
results show that the observed small-world property is 
mostly attributable to increased clustering with respect 
to corresponding random graphs (note the similarity be- 
tween Figure 2 and Figure 3). In most of the studied 
parameter space, A is very close to 1 (i.e. L similar to 
that of a random graph) ; lower values appear in general 
for parameter values that lead to disconnected graphs, 
typically with more components than random graph of 
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FIG. 3. Average values of the relative clustering coefficient 
7 of the binarized functional connectivity matrix FC gener- 
ated by AR model with random structural connectivity ma- 
trix SC. In most of the parameter space, 7>1, contributing 
to small-world topology. Note the similarity to Figure 2, sug- 
gesting the critical role of increased clustering in the observed 
small- world topology of FC. Parameter values as in Figure 2. 
Logarithmic scale used both for the ordinate and the abscissa. 
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FIG. 4. Average values of the relative average path length A of 
the binarized functional connectivity matrix FC generated by 
AR model with random structural connectivity matrix SC. In 
most of the parameter space, A ~ 1, in agreement with small- 
world topology. Parameter values as in Figure 2. Logarithmic 
scale used both for the ordinate and the abscissa. See text for 
more discussion. 



corresponding density. 

Because the small-world index a is computed by com- 
parison with random networks, it is clear that the par- 
ticular choice of the random network model may affect 
the a values. While Erdos-Renyi model is still stan- 
dard for these applications, several recent studies used 
the Maslov-Sneppen model 5 that conserves not only the 
density but also degree distribution of the graph. We 
have explored the a values observed using the Maslov- 
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Sneppen random graph model using the same framework; 
the general observation was that of weakening, but not a 
complete mitigation of the effect (see Supplemental Ma- 
terial 17 Figure 4). This suggests that degree distribution 
in functional-connectivity-based graphs is indeed differ- 
ent from that of Erdos-Renyi model, yet it does not fully 
explain the specific properties of the functional connec- 
tivity matrix, in particular the increased level of cluster- 
ing. 

With other parameters fixed, we commonly observe 
the weakest effect if the resulting FC graph has the same 
density as the generating SC. This can be understood 
since in our model only the direct structural links are 
likely to generate the strongest (and quite homogeneous) 
temporal correlations. If the number of the functional 
connectivity links is set equal to that of structural links, 
the indirect and therefore weaker links are very likely 
to be discarded through the thresholding. The practi- 
cal utilization of this fact is limited due to two points. 
Firstly, the density of underlying structural connectivity 
is typically not known in practice. Secondly, the strength 
of the structural links is commonly not homogeneous. 
Then, some indirect links can become stronger than some 
weaker direct links, causing increase in the minimum a 
value. This effect was observed in separate runs of the 
simulations that used Gaussian or uniform distribution 
for the strengths of structural links (see Supplemental 
Material 17 , Figure 5 and 6). 



IV. DISCUSSION 

We have shown that the approach of constructing net- 
work connectivity graphs from correlations of activity 
time series of local nodes leads to particular graph struc- 
tures, characterized by increased clustering compared to 
common random graph models. This may lead to at- 
tribution of small-world properties to networks possess- 
ing purely random structural connectivity. By extensive 
mapping of strength of the effect as function of model 
parameters, we documented that this phenomenon is not 
a special theoretical case or a negligible effect, but is a 
rather pronounced and general phenomenon. While the 
coverage of parameter space was necessarily limited, it is 
sufficient not only to demonstrate the existence of the ef- 
fect, but also its potential strength and existence of com- 
plex modulation by system parameters. The strength 
of the effect suggests it is of a scale relevant for inter- 
pretation of real-world results, considering that reports 
of small-world properties of real-world networks show a 
values comparable to those observed in our model - from 
the order of several units in smaller, relatively dense net- 
works (e.g. a ~ 1.5 in 18 ) to several tens in larger and 
less dense networks (see e.g. 19 ). Note that reports of 
small-world network topology based on functional con- 
nectivity analysis of measured time series is becoming 
increasingly used not only in brain imaging (see 20 and 
references therein), but also in other fields such as cli- 



mate science 3 ' 21 or economics 22 . 

As one of the most widely discussed concepts in the 
emerging field of complex network analysis, the small- 
world phenomenon has already received some attention 
regarding the methodology of its quantification. The 
two most notable attempts in this direction focused on 
the confounding effects of inadequate spatial sampling 6 
and finite-size temporal window estimation of correlation 
coefficient 7 . Some other reasons for a careful interpreta- 
tion of graph-theoretical analysis of functional connectiv- 
ity matrices have been noted in recent reviews 15,23 ' 24 . 

In the current study we bring a direct theoretical 
argument supported by extensive quantitative simula- 
tions proving that the interpretation problems regard- 
ing correlation-based small- world networks stem from the 
intrinsic properties of commonly used functional connec- 
tivity measures, rather than merely from inappropriate 
sampling. Even without the finite sample size problem, 
and with no spatial oversampling, functional connectiv- 
ity matrices of randomly coupled trivial linear processes 
show small-world properties. 

The current study utilizes linear statistical dependence 
measure (namely Pearson's correlation coefficient) and 
linear autoregressive model to study the consequences 
of the functional connectivity approach to construction 
of network graphs. A list of parameters of the model 
was varied, covering a large portion of the parameter 
space. Nevertheless, the computational requirements mo- 
tivated some limitations of the parameter swipe. The in- 
creased clustering was documented across a multidimen- 
sional range of parameter values and for several choices of 
random distributions of structural links. While the effect 
was ubiquitous, its strength depended substantially on 
parameter values - in particular, the highest values were 
typically obtained for sparse networks, an observation in 
line with results reported in literature (see 25 ). The accu- 
rate theoretical description of the parameter dependence 
of the strength of the effect for a reasonably general class 
of systems is a subject of future work, as it might help 
taking the effect properly into account in data-analysis. 

Indeed, the AR(1) process analyzed here may serve 
only as a simplest first approximation for some real- world 
systems studied with the above discussed functional con- 
nectivity approach. We conjecture that since the origin 
of the phenomenon is in the partial transitivity prop- 
erty of correlation matrix of a dependent variable set, the 
details of the stochastic process and even the fact that 
the data are generated by a temporal stochastic process 
arc not central to the phenomenon. In particular, the 
potential for spurious small-world observations may be 
also relevant to correlation-based networks constructed 
from other than time series data. Moreover, based on 
the elementary theoretical argument we conjecture that 
the main qualitative result generalizes to other depen- 
dence measures — to the extent to which these alterna- 
tive measures hold a similar partial transitivity property 
as the linear correlation. Moreover, if the linear corre- 
lation is applied on nonlinear systems, the phenomenon 
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should also appear — linear correlation may not be able 
to capture the dependences fully, but it will still show its 
specific properties. The detailed treatment of the related 
issues is a topic of further work. 

Importantly, in the current paper, we do not challenge 
the specific findings of small- world structure of the struc- 
tural networks underlying systems such as human brain. 
Rather, we only suggest that observations of small- world 
observations in functional connectivity, although poten- 
tially appealing, may not per se say much interesting 
about the underlying system, as these are probably ubiq- 
uitous in completely randomly connected dynamical sys- 
tems, in a manner illustrated on a simplistic example 
of a linear dynamical system in this paper. More gen- 
erally, this small-world bias is a specific example of the 
general notion that graph-theoretical properties of func- 
tional connectivity matrices should be interpreted differ- 
ently than those of structural matrices. 

Considering our results, in order to properly support 
the case that in a particular experimental dataset the 
observed small-world functional connectivity properties 
are significantly elevated, one would have to compare the 
experimental data with appropriate specific null-model 
(rather than that of Erdos-Renyi or Maslov-Sneppen) or 
provide a correction for the bias related to the functional 
connectivity approach. 

Unfortunately, the null model of independent linear 
stochastic processes preserving length, frequency content, 
and amplitude distribution of the original time series that 
was proposed in Ref. 7 , may not be suitable for testing 
hypothesis regarding network topology, since its alterna- 
tive includes any set of dependent processes, irrespective 
of the topology of the pattern of their dependence. Re- 
jection of the null hypothesis according to Bialonski et 
al. 7 gives statistical support for any dependence struc- 
ture from the small- world topology on one side, to the 
randomly connected processes on the other side. 

As clearly demonstrated in the simulations, the 
strength of the bias crucially depends on many param- 
eters of the underlying system, including the density of 
links and the distribution of their values. 

On the other side, researchers commonly turn to com- 
puting "at least" functional connectivity because they 
consider the full structure of the underlying dynamical 
system difficult (or impossible) to estimate from avail- 
able data. Theferore the above named parameters of the 
underlying system, necessary for determining the size of 
the small-world bias, are not available in practical situa- 
tions. This is why constructing a suitable test or correc- 
tion procedure is not straightforward and constitutes on 
open problem for further research. 

Of course, in many situations, the functional connec- 
tivity approach may still be suitable and remain in wide 
use; it is just that the interpretation should take properly 
into account the limitations of the approach. 



V. CONCLUSIONS 

In conclusion, the increasingly used functional- 
connectivity approach to complex network graph con- 
struction bears an inherent strong bias towards small- 
world-like properties as measured by the conventional 
criteria. This should be taken into account when analyz- 
ing data using this approach in conjunction with graph- 
theoretical analysis in order to avoid erroneous interpre- 
tations. Ultimately, principled solutions should be sought 
to overcome this bias as well as other complications con- 
cerning the relation between structural and functional 
connectivity 26 . Dynamical models of complex networks 
are likely to play a key role in this challenge. 
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